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Abstract 

Angular distribution measurements of ^H(^Be/Be)^H and ^H(^Be,^B)n reactions at 
Ec.m. ~ 4.5 MeV were performed to extract the astrophysical S'i7(0) factor using the asymp- 
totic normahzation coefficient (ANC) method. For this purpose a pure, low emittance ''Be beam 
was separated from the primary ^Li beam by a recoil mass spectrometer operated in a novel mode. 
A beam stopper at 0° allowed the use of a higher ^Be beam intensity. Measurement of the elastic 
scattering in the entrance channel using kinematic coincidence, facilitated the determination of the 
optical model parameters needed for the analysis of the transfer data. The present measurement 
significantly reduces errors in the extracted ^Be(p,7) cross section using the ANC method. We get 
Si7 (0) = 20.7 ± 2.4 eV b. 

PACS numbers: 25.60.Je, 25.60.Bx, 26.20.+F, 26.65.+t 
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Recently, the Sudbury Neutrino Observatory (SNO) group has discovered, in a simulta- 
neous measurement of the electron neutrino flux and the sum of three active neutrino fluxes, 
that a large fraction of the high energy electron neutrinos emitted in the /3~^ decay of 
in the Sun is transformed into other active neutrino flavors on their way to detectors on 
the Earth m. Together with the observation that the reactor produced antineutrinos also 
oscillate it seems that a solution of the solar neutrino problem is at hand. To ascertain 
more accurately if there is a transformation of solar electron neutrinos into sterile neutrinos 
the precision of both experimental measurements and theoretical predictions must be im- 
proved j^. In order to make more accurate theoretical predictions of the neutrino flux, 
the rate [or the related zero energy astrophysical 5'-factor, 6*17(0)] of the reaction ^Be(p, 7)^B 
that produces ^B in the Sun must be better determined jj]. 

Recent precision '^Be{p, 7) measurements have yielded values of <S'i7(0) which are clustered 
around 18.5 eV b P| and 22.0 eV b and are not consistent with each other within the 
quoted errors. In view of this discrepancy the determination of 5*17(0) by other methods, with 
different systematic errors, is necessary. One kind of indirect measurement uses the Coulomb 
dissociation of ^B where the latest experiment has yielded a 6*17(0) of 18.6 eV b Qj. 
The other technique uses the (''Be,^B) transfer reaction to extract the magnitude of the 

asymptotic radial wave function, characterized by the asymptotic normalization coefficient 

n 

(ANC), of the proton in ®B. This is then used to calculate the value of 5*17(0) 9]. The ANC 
method has been experimentally validated through the agreement of the measured ^^0(p,7) 
cross section and that derived from the data on ^^0('^He,d)^''F reaction [m]. Recently, the 
5i7(0) has been extracted from the transfer reactions ^°B(^Be,^Bj^Be and ^^N('^Be,^B)^^C 



yielding values of 17.8 ± 2.8 and 16.6 ± 1.9 eV b, respectively 

Because of the precise knowledge of the proton wave function in the deuteron and a 
simpler exit channel (which makes it more amenable to continuum discretized coupled chan- 
nel (CDCC) methods Q|), the d(^Be,^B)n reaction may still be an attractive choice for the 
extraction of 5i7(0) with the ANC method provided the data are taken at beam energies 
where this reaction is peripheral. It has been shown 1^ that for entrance channel center of 
mass (crn.) energies around or below 6.0 MeV this condition is fulfilled for this reaction. 



Liu et al. 



14| have reported a 5i7(0) of 27.4 ± 4.4 eV b from the analysis of t 



experiment done at i?c.m. = 5.8 MeV. However, it has been pointed out 



leir (i(''Be,^B)n 
13, Q that there 



may be a large uncertainty in the value of 5i7(0) extracted from this experiment due to lack 
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of the knowledge of the entrance and exit channel optical model parameters (OMP) as the 
corresponding elastic scattering data are not available. 

The present work is a significant improvement over the earlier work of Liu et al. 14 1. 
The ^Be beam was produced by operating the existing recoil mass spectrometer, HIRA 16| , 
in a novel optical mode leading to a beam of superior quality (purity > 99.9%, beam spot 
size ~ 3 mm, angular divergence of ±1°). This allowed the use of a stopper in the beam 
path enabling the use of a higher intensity (by a factor of ~ 10) beam thus reducing the 
statistical errors significantly. In a separate experiment, to obtain the entrance channel 
OMP, the elastic scattering angular distribution was measured using kinematic coincidence 
to eliminate the background from other target elements. These have allowed a more precise 
measurement of the d(^Be,®B)n angular distribution at the lowest energy of i?c.m. = 4.5 MeV 
and the determination of <5'i7(0) by the ANC method. 

The measurements were made at Nuclear Science Centre, New Delhi, using a radioactive 
ion beam of 21.0 ± 0.5 MeV ''Be incident on a 1.05 mg/cm^ deuterated polyethylene (CD2)n 
target. The ^Be was produced through the p(^Li,^Be)n reaction at E(^Li) = 25 MeV using 
a pulsed ^Li beam (intensity ~ 3xl0^°/sec) from the 15UD Pelletron, at a 4 MHz repetition 
rate and FWHM ~2 nsec, to obtain an average ''Be intensity of 3000/sec. The deterioration 
of the production target, a 20 fim polyethylene (CH2)n foil, was minimized by automated 
linear and rotary motions. The forward going ''Be particles were selected using the recoil 
mass spectrometer, HIRA, operated in a new ion optical mode optimized for such inverse 
kinematic reactions jjj]. In this mode the reaction products of interest were focused through 
a slit placed at the center of the magnetic dipole to reject the primary beam. The selected 
secondary ion beam was refocused at the target such that the beam spot was a replica of 
that at the production target. Two silicon telescopes placed at ± 30° with respect to the 
beam direction were used in the primary target chamber to measure the recoil protons. The 
ratio of counts in the recoil proton peak from the production target to the counts of ^Be 
in the secondary reaction chamber was monitored and kept constant to ensure an identical 
and reproducible trajectory of ''Be through HIRA. A 3 mm diameter graphite collimator 
placed 86 mm upstream of the target was used to limit any unforeseen wandering of the 
beam spot. The X-Y profile of the ''Be beam was monitored using a multi-wire proportional 
counter (MWPC) placed ~ 90 mm behind the target. The MWPC had entrance and exit 
windows of 1.5 /im polyethylene and was operated using isobutane gas at a pressure of 
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5 mbar. A 4 mm diameter tantalum disk mounted on a thin (0.25 mm 0) wire, at a distance 
of 114 mm from the target, stopped the main beam. This reduced the beam flux incident 
on the downstream detector telescope by a factor of ~8. The detector telescope consisted 
of a AE gas ionization chamber (IC) followed by a 50x50 mm^ two dimensional position 
sensitive Si-E detector (PSSD). The IC had an opening of 55x55 mm^, active depth of 
60 mm, and was operated at a pressure of 50 mbar. The position and energy resolutions 
of the PSSD and IC were measured to be <2 mm, 200 keV and 200 keV, respectively, 
using a ^^^Am alpha source. The energy, position, pileup parameter (generated from the 
zero crossover time of the bipolar E pulse) and time of flight (TOP) with respect to the 
timing reference of the beam pulsing system together with the scaled down monitor detector 
energy signals were recorded in an event by event mode. The MWPC energy output was 
recorded independently in a multichannel analyzer to minimize the dead time in the data 
acquisition system. This was used to obtain the integrated ^Be beam intensity. A 10 Hz 
precision pulser was used to monitor the dead time of the data acquisition system as well 
as the electronic gain of, and noise in, the detector system. The response of the detector 
telescope was continuously monitored using a weak ^^^Th alpha source. The count rate in 
the detector was kept constant to within ± 10% so as to keep the pileup fraction similar for 
different runs. 

The beam profile was maintained within ± 0.25 mm during the runs for an accurate 
angular definition. The stopper allowed a ten fold increase in the ^Be intensity as compared 
to that of Liu et al. for a similar pileup rate. HIRA was rotated to 2° and the scat- 
tered ^Li, ^Be and ^^C ions were selected to calibrate the particle identification (PID) of the 
detector telescope in situ. This was essential in choosing the AE-E two dimensional (2D) 
gates for ^B in conjunction with a SRIM calculation. The (CD2)„ and (CH2)„ data 
with (without) stopper were taken for 4x10* (8x10'') and 1.5x10* (2x 10'') '^Be incident 
particles, respectively. The PID parameter was calculated in the standard way using the 
AE-E information [l^. Typical PID spectra after gating on the pileup, time-of-fiight and 
particle energy (to remove the elastically scattered ^Be from the target) are shown in Fig. 1. 
The *B events were selected using suitable gates on PID, pileup parameter and TOP. The 
angular distribution was obtained using the position information from the PSSD. The po- 
sition response of the latter was measured offline using an alpha source and a mask placed 
on the detector. An overall angular resolution of 0.9° was estimated taking into account the 
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angular divergence, transverse beam profile, angular straggling in the target and gas, and 
the position resolution of the PSSD. The yields were corrected for small dead time losses 
and transmission loss in the MWPC 

Elastic scattering for the ''Be + d system was also measured at 20.3±0.5 MeV using 
the same (CD2)„ target. The schematic of the experimental setup is shown in Fig. 2a. In 
an elastic scattering event, the recoiling deuteron is detected in the annular detector (A2) 
and ''Be in the forward detector telescope AE(gas)-Si 2D. The ''Be particles reaching the 
detector telescope as a result of scattering from the collimator and ^^C (in (CD2)„ target) 
and from the beam halo have no coincident events in A2 and were eliminated. Although the 
inelastic contribution from the 429 keV state in ''Be could not be resolved, its contamination 
to the elastic scattering yield is expected to be small (also supported by low energy ''Li+^^C 
scattering measurements) 2^. This has been checked by performing a calculation for the 
inelastic excitation cross section to the first excited state in ^Be as discussed below. Fig. 2b 
shows the experimental elastic differential cross-sections. 

These data were fitted in a standard optical model analysis using the code 
SN00PY8Q The depths of the real and imaginary parts of the d-OMP were con- 

strained by the method described in Ref. [22] where the d-OMP was calculated by folding 
the nucleon optical potentials corresponding to half the deuteron energy so as to resolve the 
discrete ambiguity in the d-OMP. Following Ref. [2^, the radius parameters of both real 
and imaginary parts were varied in the range of 3.5 to 4.5 fm. In the fitting procedure, 
we applied an additional constraint that the total reaction cross section (ctr) be close to 

n 

~ 1.2 b which is obtained by calculations done with the OMP of Ref. j2^ for the rf+^Li sys- 
tem at comparable energies and also from the measurements of reported for the d+^Be 



system 



2J|. Four sets of best fit potentials obtained from a minimization analysis of 



the data are shown in Table I. The fit to the elastic angular distributions obtained with 
potential SI is shown in Fig. 2b (dashed line). Also shown in this figure is the sum of the 
elastic and inelastic (to the first excited state in ''Be at 427 keV) cross sections (solid line). 
The latter has been calculated using the same set of d-^Be OMP and p2 of 0.6 j2^. Results 
obtained with sets S2 - S4 are similar and cannot be distinguished from these curves. This 
figure also shows that the potentials sets 1 and 2 of Ref. provide very poor fits to our 
elastic data. Similar poor fits are obtained using the potential sets given in Refs. jisl . Iisj ] . 
While the fits could be improved upon by performing measurements with better statistics 
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and having more data points, the present data on elastic scattering angular distribution can 
clearly discriminate between the different deuteron optical potentials and the potential sets 
extracted by us are the only ones which provide any reasonable fit to these data. 

The measured d(^Be,^B)n angular distribution (shown in Fig. 3) has been analyzed 
within the finite range distorted wave Born approximation (FRDWBA) using the code 
DWUCK5 Q (with full transition operator) which has been modified to include exter- 
nal form factors (FF). The FF for the deuteron-neutron overlap has been obtained from 
the deuteron wave function (including both s and d states) corresponding to the Reid soft 
core potential. A two body model for the p - ^Be system has been assumed where the 
proton occupies a single particle state n£j and the ^B {p - ''Be) overlap function is written 
as S^^'^Unij{r). Here Unij{r) is the normalized single particle radial wave function and S is 
the spectroscopic factor which is directly related to the ANC P and subsequently to the 
iS'i7(0) factor. While the results obtained with the proton in the 0^3/2 orbital are given here, 
thosecalculated using the Op„, proton conBguration separately are almost identical. Five 
sets 1271 of the neutron OMPs were used in the FRDWBA transfer calculations. These were 

n 

obtained from the global parameterizations given in Ref. |2S| (used extensively in Hauser- 
Feshbach calculations), from fits to n + ^°'^^B scattering at 9.72 MeV (two sets) and to 
p +^ Be scattering at 5 and 6 MeV (two sets). The compound nuclear (CN) contributions 
were calculated using a Hauser-Feshbach code HAFEST 129^. As can be seen from Fig. 3, 
these contributions are small (~ 7.5% at 44° and ~ 1% at 8°). The uncertainty in the 
CN contributions was estimated to be around 50% (which, however, adds only about 1% 
to the overall theoretical uncertainty). This has been included in the systematic error of 
the extracted 5'i7(0). The calculated transfer angular distributions have been folded using 
a Monte Carlo simulation which took into account the spatial, angular and energy spread 
of the beam, the finite thickness of the target and the position resolution of the detector. 
The measured transfer angular distribution below 45° was used to extract the iS'i7(0). This 
was done to (a) minimize the error arising from the contributions of CN and higher order 
processes which may affect the extracted 5'i7(0) and (b) use the forward angle data to the 
maximum extent in order to reduce the statistical errors. These theoretical calculations at 
forward angles {9c.m. < 45°) were scaled to find the best fit to the data (from which the 
CN contribution is subtracted) yielding S which has been used to calculate >S'i7(0) by the 
procedure discussed in Ref y. While calculating the iS-factor the correction arising from 
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the p—^Be scattering length has been included Ssj ]. 

We have verified that the peripheral condition of the transfer reaction is fulfilled at our 



beam energy by two ways: (1) by using 
the p-'^Be system (given in Refs. jsy, |3] 



bur different sets of bound state potentials for 



32 



331 ]) which lead to a variation of the ANC 



by only ±3.5% whereas the calculated transfer cross sections changed by about 30%, and 
(2) by introducing a lower cutoff in the radial integrals where it was found that results (at 
the forward angles) obtained with a lower cutoff of up to 4.0 fm were almost identical to 
those obtained with no lower cutoff. It should also be noted that the multi step processes 
(inelastic excitation + transfer, breakup fusion) and core excitation in ®B have negligible 
effect on these cross sections |34| . 

The value of 5*17(0) and the systematic errors arising from the uncertainties in the d+'^Be 
OMP, n+^B OMP parameters and the bound state wave function of the proton in have 
been estimated from 4, 5 and 4 choices, respectively, for each of these inputs. Each of 
these 80 combinations was used to derive the best fit to the experimental transfer angular 
distribution and hence the corresponding 5*17(0). The range in which our calculated angular 
distributions lie can be seen from Fig. 3 (short and long dashed lines). The calculated 
angular distributions agree well with the data within the experimental uncertainty. The 
inset of Fig. 3 shows a histogram for the number of occurrences/0.5 eV b (N) for the 5*17(0) 
which ranges from 18.8 eV b to 22.1 eV b with a calculated mean and RMS deviation of 
20.7 eV b and 0.9 eV b, respectively. The distribution of these derived 5-factors should 
give a reasonable estimate of the theoretical uncertainties considering the large number 
of combinations for the potentials used. The total systematic error after including the 
uncertainty in target thickness (±2% by weighing samples from the same stock of target 
material) is ± 1.0 eV b. The statistical error estimated from these fits is ± 1.4 eV b. 

If the above exercise is carried out using the first 3, 5 and all 11 data points, starting from 
the most forward angles, the extracted 5*17(0) turns out to be 22. 3± 2.7 eV b, 22. 7± 1.9 eV b 
and 18. 5± 1.6 eV b, respectively. The mean S-factor from the analysis using these data 
sets and the 8 data point set used earlier is 20.6 eV b. We may use this to estimate an 
additional error arising from the different choices of data points. This is probably a very 
conservative estimate of the error since it is expected that using the larger angle data makes 
it prone to contributions from higher order processes and uncertainties in compound nuclear 
contributions while use of only the most forward data points increases the statistical error 
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while not making optimal use of the data. Nevertheless if this spread of 1.7 eV b is added 
in quadrature to the statistical and systematic errors mentioned earlier the overall error 
increases to 2.4 eV b. Since the systematic error alone contributes about 70% to the total 
error there is scope for reducing it. This would require a higher statistics elastic scattering 
and transfer measurement covering a larger angular range and the more elaborate CDCC 
calculations. 

In conclusion, we have measured, for the first time, both the (i(^Be,^B)n transfer and 
the entrance channel elastic differential cross sections at the lowest beam energy hitherto 
using a high quality ''Be beam from a recoil mass spectrometer operated in a novel optical 
mode. The extracted ^Be(p,7) S'i7(0)-factor is 20.7 ± 2.4 eV b is in good agreement with 
the latest direct capture measurements |S| and those determined from the CDCC analysis 
of the ^B breakup reaction ^3^. Thus the disagreement in the values of 817(0) determined 
by direct and indirect methods is reduced (see also 37|). This experiment, therefore, clearly 
demonstrates that the ANC method can be used for reasonably precise measurements of 
other (p,7) ^-factors involving short lived nuclei, where direct capture measurements may 
be very difficult, if not impossible. 
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TABLE I: Parameters of the Woods-Saxon optical model potentials extracted from the analysis 
of the present d + ^Be {Ec.m. = 4.4 MeV) elastic scattering data. A spin orbit term with Vgo = 
8.60 MeV, rso = 2.17 fm, and ttso — 0.6 1 fm, has been added to all the potential sets. The optical 
potential is defined as that in Ref. [^^I with the light convention for the radius. 



Pot. 
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SI 


103.12 


2.23 


0.62 


79.03 


2.37 


0.17 


S2 


107.87 


2.17 


0.61 


58.84 


2.28 


0.25 


S3 


92.54 


2.41 


0.57 


117.50 


2.45 


0.14 


S4 


121.49 


1.97 


0.66 


54.88 


2.38 


0.28 
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FIG. 1: a) Particle identification (PID) spectrum for the polyethylene target with suitable energy, 
pileup and TOP cuts (see text) for a integrated ^Be flux of 1.65x10^ particles, (b) PID spectrum 
for deuterated polyethylene target. 
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FIG. 2: a) Schematic of detector setup used for elastic scattering measurement, b) Angular 
distribution for the elastic scattering of deuteron on ^Be at Ec.m. = 4.4 MeV. Curves are fits to the 
data using optical potentials SI (solid), and S4(dotted) of Table I. Results obtained from potential 
sets of S2 and S3 are not distinguishable from those shown. The dashed and dot-dashed lines 
represent the cross sections obtained with potential sets 1 and 2, respectively of Ref. [l^ . 
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10 20 30 40 50 60 

FIG. 3: Measured d(^Be,^ B)n angular distribution together with the folded FRDWBA + CN 
cross sections shown as solid and short and long dashed lines lines (see text). The calculated 
compound nuclear contributions are shown by the dotted line. The inset shows a histogram plot 
of the extracted 817(0) using the various combinations of d-^Be, n-^B OMP and p-^Be bound state 
potentials (see text). 
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